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On a Graphical Method for determining the Orbit of a Binary Star. 

By Professor S. Grlasenapp. 

(Communicated by Herbert Sadler .) 

The problem of determining the elements of a binary star 
consists of two separate parts : (1) to obtain the apparent orbit 
from a series of given observations, and (2) to compute the ele¬ 
ments of the true orbit, the apparent orbit being given. Both 
parts can be solved by graphical and analytical methods. 

Between the true and apparent orbit exists a geometrical 
relation which furnishes the means for obtaining one orbit when 
the other is given; this is a purely geometrical problem. 

The apparent orbit is never known; we can only determine 
the probable apparent orbit of the component stars. Each com¬ 
puter tries to obtain the most probable orbit. 

The general equation of an ellipse, being a curve of the second 
degree, is— 

ax + fry + 7a’ 2 + Bxy + ey 2 4-1=0.(1) 

where x and y are the rectangular co-ordinates of the comes 
relative to the principal star which coincides with the origin of 
the co-ordinates. The plane of co-ordinates is perpendicular to 
the line of sight. 

The coefficients of the equation (1) for an ellipse must satisfy 
the well-known conditions— 

(fid — 2ae) 2 — ( 8 2 — 4 €7) (fi 2 — 4c) > O, 

.(2) 

5 2 — 457 <0. 

I should add to them the two following :— 

7 < o and e < o,.(3) 

which represents that the principal star lies in the inner part of 
the ellipse—i.e. that the ellipse obtained intersects the two co-or¬ 
dinate axes of x and y on opposite sides of the principal star. If 
the conditions (3) are not fulfilled all the parts of the ellipse 
lie outside the principal star. In this case, in determining the 
elements of the true orbit by the usual methods, we shall obtain 
imaginary elements. 

Let us now consider a very simple method for determining 
the coefficients a, /3, y, 8, and e of the general equation (1). 

Let us assume that the apparent ellipse is drawn with such 
care as is generally employed in graphical methods. 

If we put y— o. the equation (1) will be reduced to an equa¬ 
tion of the second degree in x — 

ax + yx 2 +1=0.(4) 

whose roots are the co-ordinates of the two points where the 
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apparent ellipse intersects the axis of x. These roots being x 1 
und x 2 , we obtain 

.•••••• ( 5 ) 


Xi +Xo I 

«= -- 1 — 7 = — 


DO-^OOty Vb-^* 27*7 

If, further, we put $=0, we obtain an analogous equation 
in y— 

% + e/r 4 - 1 = 0 , ..(6) 


whose roots are the co-ordinates of the two points where the 
apparent ellipse intersects the axis of y. These roots being y Y 
and y 2 , we obtain 


is — 

ViU-i 


I * 


( 7 ) 


If, therefore, we measure with compasses the co-ordinates aq, 
a? 2 , 2/1, and the four intersecting points of the apparent 

ellipse with the co-ordinate axes, we obtain immediately, by 
means of the formulae (5) and (7), the values of the four coeffi¬ 
cients a, / 3 , 7, and e. 

For the determination of 8 we must take a fifth point. Let 
its co-ordinates be x s and y 3 ; then, from the general equation (1), 
we obtain 

§ _ __ 1 + ax s + jjjga + 7*s 2 + nhf ....... (8) 


From the form of this equation we can see that the most 
accurate determination of 8 will be obtained in the case when the 
divisor x 3 y 3 has its maximum value. It is always possible to 
choose, without any calculation, such a point on the apparent 
ellipse whose co-ordinates very nearly satisfy this condition. As 
a check we can take a second point on the opposite side of the 
apparent ellipse, and obtain a second value of 8. 

The mean of both values will be more accurate and more free 
from accidental errors which may have been made in drawing 
the apparent ellipse. 

It is easy to see that the calculation of the coefficients a, f 3 , 
y, 8, and e by means of the formulse (5), (7), and (8) is very 
simple. The results are quite independent of the errors which 
may occur in drawing secondary ellipses, circles, and lines. This 
is an advantage which is worthy of attention. Further, all the 
calculation can be made in 10 to 15 minutes of time. 

To illustrate this method, I will here give a numerical 
example. I take the excellent “ Handbook of Double Stars ” of 
Edw. Crossley, F.H.A.S., Joseph Gledhill, F.R.A.S., and James 
M. Wilson, M.A., F.R.A.S., and select the example of a graphical 
determination of the orbit of Castor. 

On Plate II., p. 116, of this “ Handbook,” we find the apparent 
ellipse of Castor. From this drawing we take the following co¬ 
ordinates of the four points of intersection of the apparent 


* The co-ordinates y x and y 2 , and also x x and x. z , must have opposite sings ; 
therefore 7 and e must he negative. 


Y 2 
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ellipse with the axes of x and y. The positive end of the axis of 
x is directed to the north (0=o°), the positive end of the axis, 
of y is directed to the point whose position is determined by 
0 — 90°. 

We have— 

for y = o x l = + 11 l mm, 8 x 2 = — 233 lr * m, o, 

for x = o y 1 = + 170 mm> 4 y 2 = - 115 mm *6, 

whence we obtain— 

loga = Iog —= 7-6676?^— 10, log)S = log T—=74444—10 

L *1*2 J L V\V 2 J 

log 7 = log j—1 = 5-5842^-10, log e = log j J —1 = 57056,*-io. 

L^i^J LytfzJ 

To obtain 8 we choose a point which is determined by the 
co-ordinates— 

^ 3 = ~ 20o mm y 3 = - 75 mm , 

which gives os— 

iogs=io g r— 1 ^ ++ ^ Xs + e ^ 3 ~i = 4^230—10. 

L J 

To check this we take another point on the opposite side of 
the ellipse, whose co-ordinates are— 

# 4 = + 5o mm y 4 = + 15o mm , 

and thus obtain a second value of log 8=4*8549. The logarithm 
of the mean of these two values of 8 is equal to— 

log 8 = 4-8392—10. 


Thus we have obtained the values of the coefficients a, ft, y, 
8 , and c of the equation (1). 

This example is a good proof of the facility with which the 
coefficients a, ft, y, &c. can be obtained by the proposed method. 

In investigating the true orbit I always use the very elegant 
formulae of the late Russian Professor of the Imperial University 
at Kazan, M. Kowalsky,^ which give the following relations 
between the elements of the true orbit and the coefficients a, ft, y, 
&c., of the general equation (1):— 


tan 2 i . - s- i r, 

— . sin 2Q = 0 —3 


q 

tan 2 i 


r 


. cos 2£1 = (7 — e) — (a 2 — /3 3 ). 


, + w i= _ (7+f)+i(a2+ ^ 

e sin A = — 2 . (/3 cos Q — a sin ft) cos i. 


• ( 9 ) 


e cos A= — 1 (0 sin Q +a cos Q). 
2 


a-- 


1 — e 2 


* These formulae may be found in the Proceedings of the Kazan Imp. 
University for the year 1873. 


% 
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•where— 

i is the inclination of plane of the orbit to the line of sight (plane :-f xy). 
£> is the position-angle of the node. 

■ e is the eccentricity. 

A is the angle between the node and the periastron ( 7 r — <3). 
a is the semi-axis major. 


The first and second equations determine the node and 
—: from the third we obtain —, and hence q and i ; then the 

<f r 


fourth and fifth equations determine e and A, and from the last 
one we obtain the value of a. 

The calculation of the elements by these formulae, as may 
• be seen, is very simple. The geometrical elements can be obtained 
iln half an hour when a, / 3 , y, S, and e are given. 

For example, we deduce— 


Further— 


.and hence— 


log . sin 2 52 = 5-1264 

log . cos 2 S 3 — 4^9464 

q2 

log tan 2a = 01768 
29 > = 56°*35 

log 5-2060 
q2 



and finally— 


lo g ~ = 5 ' 9°55 
2 2 

log q = 2-1978 
log tan 2 i = 9-601 5 
log e sin A = 9-491 i n 
log e cos A = 9 3420 
log tan A = 0’i49i ft 
log 6 = 9-5901 


loga = log-i~ 2 


2-2816 


S3 =28°-i8 


i= 32 c ‘ 3 ° 


x = 3 ° 5°'35 
c = 0-380 


a = 184™"” 3 


To convert the value of a into seconds of arc, we take from 
the “ Handbook,” on p. 113, the distances on the millimetre scale. 
We find there that— 

1238“- 2 = 45-"o7 
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Therefore— 

i84 nmi, 3 = 6"7i. 

Let us compare the geometrical elements of Castor thus 
obtained with those which are found by the authors of the 
“ Handbook ” :— 

Elements of the Handbook. Glasenapp's Elements. 


a — 

6'67 


671 


0*38 

e *s» 

0-38 

Pc ■« 

28°25 

£ = 

28-18 

i — 

32-25 

i = 

32-30 

X — 

305-17 

X = 

305-35 


The elements T and n must be determined by one of the usual 
methods. 

Observatory of the Imperial University } 

St. Petersburg. 


On the Determination of Normal Places. By Lieut-General 
J. F. Tennant, R.E., F.R.S. 

When it is proposed to correct an approximate orbit by 
comparing the places deduced from it with a large number of 
observations, and obtain corrections to the elements by the 
method of least squares, each observated element of the place 
would produce an equation of condition. It becomes then 
necessary to reduce the number of these equations of condition 
somehow, so as to limit the amount of computation in forming 
normal equations. 

It seems to have been the practice to compute the equation 
of condition for each observation, and then to take the mean of 
the equations corresponding to a group of observations, and use 
it as one equation. How the use of normal places formed by 
comparing the observations themselves with an ephemeris from 
approximate elements is practically universal. Of course the 
accuracy of these normals is vital. Ho care in computing can 
produce trustworthy results unless the normals are accurate, and 
to make them so a considerable number of observations have to 
be used, which probably extend over a considerable time. The 
determination of normal places has been treated of in the fourth 
volume of Oppolzer’s Lehrbuch , but his remarks do not seem to 
me to quite meet the wants of those who have not all his experi¬ 
ence and skill to guide their steps. 

The difference between the series of places computed from 
the approximate elements and the true ones may be expressed hy 
a series of the form 

a + bt 2 + ct~ + &c. 
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